Do gluons carry half of the nucleon momentum? 
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We examine the conventional picture that gluons carry about half of the nucleon momentum in 
the asymptotic limit. We show that this large fraction is due to an unsuitable definition of the 
gluon momentum in an interacting theory. If defined in a gauge- invariant and consistent way, the 
asymptotic gluon momentum fraction is computed to be only about one fifth. This result suggests 
that the asymptotic limit of the nucleon spin structure should also be reexamined. A possible 
experimental test of our finding is discussed in terms of novel parton distribution functions. 
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A classic prediction of perturbative quantum chromo- 
dynamics (QCD) is that about half of the nucleon mo- 
mentum is carried by gluons in the asymptotic limit [H- 
This renowned fraction has become a deeply rooted pic- 
ture in the minds of hadron physicists, partially because 
it was supported 0, indirectly, by measurement of the 
quark momentum fraction in deeply inelastic scatter- 
ing (DIS) of leptons off nucleons (which we will com- 
ment on later). Considered together with the nucleon 
spin problem discovered later, however, the taken-for- 
granted knowledge about the gluon contribution to the 
nucleon momentum becomes questionable: Were the 
gluon momentum known perfectly, why was severe dif- 
ficulty (namely, with gauge invariance) encountered with 
even a sound theoretical definition of the gluon orbital 
angular momentum, which should be a direct extension 
of the gluon momentum? 

To better appreciate and clarify the issue, let us recall 
that the conventional gluon momentum fraction is based 
on the following decomposition of the total momentum 
operator in QCD: 



^total — 



J d^xtP^'^Di^ + J (fxE X B = + ,^g. (1) 



Here D — V — igA is the covariant derivative. (We 
suppress all color indices and generators.) The quark 
and gluon momentum fractions are defined through the 
nucleon-state expectation values of the renormalized op- 
erators and QCD prescribes the scale evolution 

of and according to the following well-known 
mixing matrix at leading order [l[: 



^ 47r 



(2) 



which leads to the well-known asymptotic limit, = 
-Ptotah where rig = 8 is the number of gluon fields 



and Uf is the number of active quark flavors. The most 
typical case oi Uf — 5 gives ~ ^Aotai, which says 
that gluons carry about half momentum of a nucleon (or 
virtually any hadron) in the asymptotic limit. 

The doubtful point that we noted above is that the 
gluon momentum density £^g{x) = E x B (the Poynting 
vector) in Eq. ([T]) cannot be used to construct a gluon 
orbital angular momentum. In fact, in a pure Yang-Mills 
theory (namely, without quarks) J cPxx x 3^g{x) gives 
the conserved total angular momentum of the gluon field, 
including both spin and orbital contributions: 



(fxExA+ I (fxxxE'VA\ (3) 



(fxxx{ExB) 



If we insist that a theoretically sound momentum den- 
sity P{x) should match the orbital angular momentum in 
the standard form j d^xx x P{x), then Eq. ([3]) indicates 
that the gluon momentum density should be identified 
as Pg{x) = E^VA^. It has often been argued that there 
is no essential difference between Pg{x) and ^g{x) since 
they give the same integrated gluon momentum. This 
equivalence, however, is limited to the pure Yang-Mills 
case: 



Ex B = E'VA' - V'iE'A) + {V ■ E)A. 



(4) 
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The second term is a total derivative and can be dis- 
carded in integration. By the QCD equation of motion, 
the third term (where I? — V — ig[A, ] is the covari- 
ant derivative for the adjoint representation) is zero in 
the absence of quarks, but not zero otherwise! Hence we 
see that the quark-gluon interaction renders / d'^x^g{x) 
and J d^xPg{x) no longer equal. Then, naturally, we may 
find a totally different picture of the nucleon momentum 
partition if instead of / d^x3^g{x) we define / d^xPg{x) 
as the gluon momentum. This apparently interesting and 
important issue, however, has (to our knowledge) never 
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been pushed to actual calculations. The obstacle is evi- 
dent but not fully appreciated: that of gauge invariance, 
which we elaborate in some detail in the following. 

The density Pg{x) is obviously gauge-dependent. The 
gauge-invariance of its integration, J (PxPg{x), is guar- 
anteed in a pure Yang-Mills theory by the equality to the 
gauge-invariant expression j d^x3^g{x). When quarks 
are present, however, Eq. ([4]) says that / d^xPg{x) loses 
its equality to / d^x3^g{x), and thus loses simultaneously 
its property of gauge invariance. 

Being a different density, Pg [x) does not give the same 
angular momentum J d^xx x Pg{x) as J d^xx x 3^g(x) 
even in a pure Yang-Mills theory, and is always gauge- 
dependent. In fact, for several decades it was stated in 
common textbooks that the separation of spin and or- 
bital angular momentum for a gauge field is conceptually 
not possible, even for the Abelian case [Sj]. If applied 
to our momentum discussion, this seems to dictate that 
the gluon momentum in the quark-gluon interacting case 
is intrinsically not meaningful, since without a proper 
orbital angular momentum definition we would have no 
means to match a momentum expression. In a pure Yang- 
Mills theory, it is possible to ignore the gauge dependence 
of / d^xx X (i?*VA'), and define the momentum density 
as E'^VA^, which, though gauge-dependent itself, does 
give a gauge-invariant integrated momentum. However, 
as we explained in the preceding paragraph, such a semi- 
justification is lost in the quark-gluon interacting case. 

It is quite tempting to agree with the common text- 
books that in gauge theories the basic physical notions 
like the gluon or photon spin and orbital angular mo- 
mentum are prohibited by gauge invariance, because 
construction of these quantities necessarily involves the 
gauge field A. The gauge-invariance problem has long 
prohibited a consistent investigation of the nucleon spin 
structure f3|, and, fairly speaking, was simply overlooked 
in the examination of the nucleon momentum. 

Fortunately, a novel solution was achieved recently in 
Ref. which supplies a unified and gauge- invariant 
treatment of the gluon momentum and angular momen- 
tum. The key to our solution is to recognize that in the 
gauge coupling ■ip'^i^i^A'^ip, the gauge field A'' plays a dual 
role: it provides a physical coupling to the Dirac field ip, 
as well as a gauge freedom to compensate for the phase 
freedom of ip. Our idea for solving the gauge-invariance 
problem is to decompose this dual role by seeking a 
unique separation = A'^^^.^ + A'^^y^y^, with A^p^,.^ a pure- 
gauge term, transforming in the same manner as does 
the full A^^, and always giving null field strength (i.e., 

^pure — <^ ^puio ^purc ' 'i/ Impure i ^puroJ 

^phys' ^ physical term transforming in the same manner 
as does F^^" . Namely, ^phys is gauge invariant/covariant 
in Abelian/non- Abelian case, while ^(Jure has the same 
gauge freedom as A^^ and can be used instead of A^ to 

+ igAi;^ 



acting on the fundamental representation and V'^^^^ = 
— iglA'^^j.^, ] acting on the adjoint representation. The 
field separation is to express ^phys and Af^^^^ in terms of 
A^ , thus does not interfere with the gauge condition or 
canonical quantization, which manipulates the full A^. 
In this way, the spin and orbital angular momentum of 
the quark and gluon fields can all be constructed gauge- 
invariantly at the density level: 



Jtotal = / d-^xE X Aphys + / d'^XX X £;'X»pure^i 



phys 



+ I d^x^p^^^:tlj + I X ■(/'^ji5puroi/'- (5) 



Here we have improved over the expression in Ref. [j| 
by using the pure-gauge-covariant derivative for Aphys 
as well, which safely guarantees the gauge-invariance of 
-E^X'puioApjjyg. This is somewhat a natural choice, since 
in QCD ^phys is gauge covariant instead of invariant, thus 
as for ij} or F'"', a covariant derivative should be applied. 
By doing so we are left with more freedom in choosing 
a defining equation for Aphys (e.g., Ppurc • ^phys 0), 
while in Ref. Q the equation had to include [Aphys, E] — 
so as to guarantee the gauge invariance of -E'VApj^y^ 
with an ordinary derivative. We have checked that all 
consistencies in Ref. Q are maintained in the present 
formulation. Moreover, superior to [Aphys, E] = 0, in 
the perturbative region Ppm-c ■ ^phys ~ together with 
^pure ~ can give explicit series expression, in powers of 
5, for ^phys (which has trivial boundary behavior) and 
hence also A^^^^ = - Aphys- As was noted in Ref. [i], 
a proper definition of ^phys non- Abelian theory is 

quite subtle, but as we explain below, this subtlety does 
not affect the present calculation. 

From Eq. (O, we can read out the corresponding 
gauge-invariant quark and gluon momentum operators: 

/\otal = j d^Xll}'^-^Dpnrei^ + j d3a;£:»25pu„A^hys 



(6) 



construct a covariant derivative -Opm-e 



pure 



It can be proven as in Ref. [J| that Pq + Pg equals 
+ S^g in Eq. ([1]). Nevertheless, the individual terms 
are distinct in the presence of quarks: l^g differs from 
Pg by a gauge-invariant quark-gluon interaction term: 

d^xE X = y d^a;£:*25pu,cA^hys + ^ d^x^)'^ gA^x^y^^. 

(7) 

Equipped with the theoretically sound gluon momen- 
tum definition which precisely matches an equally sound 
definition of the gluon orbital angular momentum, we 
can now take up the previously untouched calculation to 
reveal a possibly different picture of the nucleon momen- 
tum partition. The calculation is straightforward and 
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parallels [l[ that based on the operators in Eq. ([T]). Let 

us first display the mixing matrix for Pg and Pgi 




which gives a different asymptotic limit for the renor- 
malized gluon momentum, P/^ = — ^% — Potai- For the 

typical case oi Uf — 5, this gives ~ l-Ptotah as com- 
pared to ~ iPtotal- 

Some technical details are worth mentioning. Because 
the total momentum is conserved, the 2x2 evolution 
matrix of the quark and gluon momenta (in whatever 
definition) has only two independent elements, as can 
be seen in Eqs. ([2]) and ([8]). The easiest way of ob- 
taining these two matrix elements is by calculating the 
expectation value of the gluon momentum operator in a 
quark state and the expectation value of the quark mo- 
mentum operator in a gluon state, which both start at 
order a^. By doing so, we avoid the task of comput- 
ing the quark and gluon wavefunction renormalization 
if we consider only the leading-order evolution matrix. 
This strategy is advantageous in obtaining Eq. ([2]), and 
becomes almost indispensable when obtaining Eq. ([8]). 
The reason is that computation of the gluon wavefunc- 
tion renormalization necessarily involves the non-Abelian 
three-gluon vertex, while for 1-loop calculation of the 
gluon matrix element in a quark state, the gluon field 
behaves like eight independent Abelian fields. Unlike in 
the non-Abelian case, the separation of an Abelian 
is unambiguous With vanishing boundary values for 
perturbative calculations, Apurc is just the longitudinal 
field y4|| = VV^^(V • A), and Aphys is just the trans- 
verse field = -V X V-2(V xA) = A- VV-2(V • A). 
(^phys ^^'^ ^puro relevant here.) Therefore the 

1-loop insertion of Pg in a quark state becomes simplest 
in Coulomb gauge W ■ A — 0, which gives Apurc = and 
Aphys = A. The calculation is then a standard textbook 
exercise, and leads to the first column in Eq. ([8]). These 
two matrix elements are j of their counterparts in Eq. 

As to the 1-loop insertion of Pg in a gluon state, we 
note that for a physical gluon, J d'^ xip^ Ap\iysip vanishes 
identically at 1-loop order. Then Eq. ^ means that 
Pg and ii^g give the same 1-loop expectation value in a 
physical gluon state, so the right-hand columns of Eqs. 
([2|) and ([8]) are the same. 

Thus we see that if the gluon momentum is defined 
properly, its contribution to the nucleon momentum is 
much smaller than conventionally accepted. The differ- 
ence can be traced to the fact that the conventionally 
defined gluon momentum actually includes a quark- 
gluon interaction term, as indicated by Eq. ((7|), while Pg 
is constructed solely in terms of the physical degrees of 
freedom of the gluon field and is thus justified as a proper 
gluon momentum. 



The novel and distinct asymptotic limit of the nucleon 
momentum partition we have presented naturally sug- 
gests that the asymptotic limit of the nucleon spin par- 
tition be reexamined according to the decomposition in 
Eq. ([5]), where the gluon spin and orbital angular mo- 
mentum are constructed (as for Pg) solely from the physi- 
cal gluon configuration. We expect that this investigation 
would also reveal a distinct picture of the nucleon spin as 
compared to that in Ref. , which was based on a light- 
cone-gauge formulation and thus included non-physical 
degrees of freedom in the gluon angular momentum. 

In the remainder of this paper, we discuss the delicate 
issue of possible experimental measurement of the proper 
quark and gluon momenta as defined by Eq. ([6]). To 
this end, it is very illuminating to first recall how the 
interaction-involving momenta in Eq. ([T]) are measured 
in hard processes. In fact, one may ask a key question: 
Why are interaction-involving quark and gluon momenta 
measured, instead of the proper ones? The answer lies, 
again, in gauge invariance. 

Collision experiments do not measure the matrix ele- 
ments of quark and gluon momentum operators directly. 
These matrix elements are extracted from the measured 
cross sections either via the operator product expansion 
(OPE) analysis of DIS, or via their relation to the quark 
and gluon parton distribution functions (PDFs) which 
can be shown to factorize in DIS and other hard pro- 
cesses Q . The reason for the conventional OPE approach 
to employ 3^g{x) and ^g{x) is evident: Expansion of the 
products of the gauge-invariant quark electric currents 
requires gauge-invariant operators. The naively proper 
momentum operators ^p^i'^ip and E'^VA^ are gauge de- 
pendent. Conventionally, the only known way of accom- 
plishing gauge invariance is by coupling a gauge field 
to ip, and by using F^"^ instead of in constructing 
physical quantities of the gauge-field. Either way, the 
physical content of the original quantity is substantially 
modified, because, as we explained earlier, a gauge field 
introduces not only a gauge freedom but also a physi- 
cal coupling. Moreover, the restriction to F^" severely 
limits one's capability in constructing certain quantities, 
e.g., the spin. Now that we have decomposed the dual 
role of the gauge field by separating its pure-gauge com- 
ponent from its physical component, we can take exactly 
the needed part and discard the unwanted part for any 
desired goal. Namely, solely to restore gauge invariance, 
we use instead of the full A^ , while Aphyg is used 

instead of F^'^ whenever the gauge-field canonical vari- 
able is unavoidable in constructing a physical quantity. 
This is exactly how we succeed in constructing the gauge- 
invariant and proper momentum, spin, and orbital angu- 
lar momentum of quarks and gluons in Eqs. fS]) and ([6|). 
In principle, the gauge- invariant Pq{x) and Pg{x) can be 
used instead of 3^q{x) and ^^g{x) in OPE analysis of 
DIS, so that proper quark and gluon momenta can be 
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measured. 

It is less straightforward to see why the conventional 
factorization approach also refers to the interaction- 
involving quark and gluon momenta. In this approach, 
what one manipulates are the quark and gluon PDFs, 
which factorize in cross sections for certain hard pro- 
cesses 0, and can be integrated to give the quark and 
gluon momenta by moment relations 0]. Naively, one 
may expect that a quark/gluon PDF should provide an 
unambiguous representation of the quark/gluon prop- 
erty, but in conventional practice this is not the case. 
The key point is how one actually constructs the quark 
and gluon PDFs. Let us first look at the quark PDF in 
a target A 

dx- 



q/A 



(0 



47r 



r 

'V exp{ig / 
Jq 



dy-A+{Q,y-,0^)}i^{Q))A, (9) 



where a path-ordered (V) gauge link (Wilson line) is in- 
serted to achieve gauge invariance. It is exactly this 
gauge link that makes the PDF defined by Eq. ([9]) an 
interaction-involving one, since the gauge field brings 
not only a gauge freedom but also a physical cou- 
pling to the quark field. The interaction term is more 
evident in the moment relation: J^oo'^^^'^i/AiO = 
2(p+)2 {i^l'^iD'^ip) A. This is just the -I- component of 

3^q{x) in Eq. dT]). Here the gauge field in £>+ originates 
exactly from the gauge link in Eq. It is often taken 
for granted that gauge invariance must be achieved at 
the price of a gauge coupling, and hence the measurable 
momenta are always the interaction-involving ones in Eq. 

But we do not really have to accept this confabulation 
for gauge invariance. Solely to guarantee gauge invari- 
ance, the gauge link can be constructed with the pure- 
gauge component A^^.^ instead of the full A^, and thus 
we can define a proper and gauge-invariant quark PDF: 



which integrates to give the desired proper quark 
momentum in Eq. ([5]): J^^,d£,^Pq/A{0 = 

Analogously, the conventional gluon PDF 



dx 



-oo '2<P^ 



(F+''(0,x-,Oi) 



xVe^ig / dy-A+{0,y-,G^)}F+{Q))A (11) 
Jo 

can be revised according to our strategy as 

P,/a{^)= r ^e-«^"--(F+^(0,x-,Ox) 



rexp{tg dy-A+{0,y-,0^)}Al^^M)A, (12) 



where, besides the pure-gauge link, the physical compo- 
nent ^phys is used instead of Fj,+ as the gauge-invariant 
canonical variable. The second moments of ^g/A and 
Pg/A relate to the interaction-involving and proper gluon 
momentum in Eqs. ^ and ([6]), respectively. 

It can be expected that the proper quark and gluon 
PDFs, Pq/A and Pg/Ai should (though probably non- 
trivially) factorize in the same processes to measure ^q/A 
and ^g/A 1 so that the proper quark and gluon momenta 
in Eq. ([6]) can again be measured! 

In closing, we remark that our approach is especially 
superior in gauge-invariant construction of polarized and 
transverse-momentum dependent PDFs with a clear par- 
ticle number interpretation, and off-forward PDFs which 
can be measured to infer the orbital angular momenta 
in Eq. (O, e.g., the polarized gluon PDF can be defined 
gauge-invariantly as 

/oo T — 

^^e-«^"^"(i^+*(o,x-,o^) 

xVexpiig [ dy-A+^^^{0,y-,0±)} 
Jq 

Xe^. + ^phy.(0))A, (13) 

with a first moment related to the gauge-invariant gluon 
spin in Eq. ([5|). Details will be reported elsewhere. 
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